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Abstract 

A global QCD analysis of the direct photon production process from both fixed 
target and collider experiments is presented. These data sets now completely cover 
the parton x range from 0.01 to 0.6, thereby providing a stringent test of perturbative 
QCD and parton distributions. Previous detailed studies of direct photons emphasized 
fixed target data. We find most data sets have a steeper pt distribution than the QCD 
prediction. Neither global fits with new parton distributions nor improved photon 
fragmentation functions can resolve this problem since the deviation occurs at different 
X values for experiments at different energies. A more likely explanation is the need 
for additional broadening of the kt of the initial state partons. The magnitude and the 
possible physical origin of this effect are investigated and discussed. 
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1 Introduction 



Direct photon production in hadron collisions is an important process to study in pertur- 
bative Quantum Chromo dynamics (QCD) due to the clean measurement of photons, in 
contrast to jet physics with the jet definition uncertainties. In addition, sensitivity to the 
initial state gluon in the important Compton scattering mechanism provides a valuable op- 
portunity to directly measure the gluon content inside the proton. These features have been 
utilized in studies ^ of earlier direct photon experiments in the context of next-to-leading- 
order (NLO) (q!Q!s) QCD calculations of Aurenche et al. 0, and in the global analyses 
of parton distributions.]^, |, The potential usefulness of the direct photon experiments 
demonstrated in these works has yet to be fully realized due to the limited range of the fixed 
target experiments as well as remaining theoretical uncertainties related to scale-dependence 
and fragmentation processes. 

In hadron colliders, earlier data from the ISR and SppS colliders have been recently 
augmented by high statistics data from the Fermilab TEVATRON. The x-range covered by 
these collider and fixed-target experiments now spans the entire interval (0.01 < x < 0.6). 
The combined range and improved accuracy of these data, as well as recent theoretical 
developments, provide a good opportunity to make a global comparison between theory 
and experiment. We report here such a study. The next two sections examine existing 
work on fixed-target and collider data and treat them in a unified context. A pattern of 
deviation of the observed Xt distribution from existing QCD expectations is highlighted. 
This is followed by a new global analysis of all existing data. We show that it is difficult to 
account for the observed pattern by new sets of parton (mostly gluon) distributions because 
the deviations occur at different x values for different experiments. Various alternative origins 
of the observed pattern and possible related phenomena in other processes are discussed. 

2 Fixed Target Results 

Previous parton distribution analyses incorporating direct photon data have focused on fixed 
target experiments, particularly the WA70 experiment (pp collisions at y/s = 23 GeV 0) 
from CERN. The ABFOW analysis performed fits on DIS data from the BCDMS exper- 
iment and used WA70 data as constraints to determine the gluon distribution, then checked 
many other direct photon experiments with the resulting parton distributions. The WA70 
experimental normalization was fixed at unity.[] The MRS analyses P, Q used a similar 
procedure but included more DIS data and added Drell-Yan results as constraints. Both 
of these groups used the "optimized scale" scheme in choosing the renormalization and 
factorization scales. Fig. [1| shows the comparison of WA70 data with the ABFOW fit in the 
form (Data - NLO QCD) /NLO QCD as a function of photon xt {2pt/\/s). The production 
is fairly central so that photon Xt is approximately the parton x value probed. One sees that 
the fit is satisfactory. To see the scale-dependence of the NLO theory, one can compare the 
results obtained from two choices of scale, for example optimized and The latter is 

normalized to the former (as are the data points) and shown in the plot as the dashed line. 

^Since the Compton cross-section is proportional to as(/i, A) • G'(x,/i), the normahzation factor directly 
affects the determination of A and G(x,/i). 
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The optimization procedure results in higher cross-sections since it yields scales smaller than 
for this kinematic region. 
In the CTEQ2 global analysis^, in addition to WA70, data from the fixed target ex- 
periments UA6 {pp collisions at y/s = 24 GeV [§) and E706 (pBe collisions^ at ^/s = 30.7 
GeV were used. These data were treated on the same footing as the DIS and Drell-Yan 
data in an overall global fit. The uncertainties on experimental normalization and theoretical 
scale-dependence are both taken into account in the analysis procedure (see Sec.| and Ref. [^] 
for details). In Fig. |^ data from these three experiments are compared to the CTEQ2M fit. 



In this figure the more recent data from UA6 (pp and pp collisions [11|), which agree with 



their earlier data, are shown. The absolute normalizations of all data and theory are used 
in this figure. These fixed-target experiments cover similar xt ranges (with E706 extending 
somewhat lower) and are roughly consistent with each other. We note that the E706 data 
show a deviation from NLO QCD on this comparison plot which was not evident in previous 
comparisons with WA70 data alone. This pattern is also seen with other current parton 
distribution sets such as MRSDO [H. 



3 Collider Results and Calculations 



The first direct photon data from hadron colliders came from the ISR pp collider ||12[, followed 



by the SppS collider |13|. These data were found to be in qualitative agreement with NLO 
theory in Ref.[|l|. However, the 1989 data from the TEVATRON pp collider experiment 
CDF [|T^, extended into a new x region (0.01), and showed a clear deviation (1.5 — 2 times 
higher) from the NLO QCD predictions below pt ~ 20 GeV. Improved parton distribution 
sets such as CTEQl p and MRSDO 0], using newer small-x deep inelastic scattering data, 
removed approximately one-half of this discrepancy, but the remaining difference persisted. 
This effect has recently been confirmed with much higher statistics data from CDF |15|. The 
obvious source of uncertainty due to choice of scale cannot be responsible for the discrepancy 
since it produces small normalization shifts with no change in shape [Q. In addition, the 
optimized scales scheme, favored in fixed-target fits, was shown to make the disagreement 
between data and NLO QCD worse |T^. Thus, much attention has been directed to the 
proper theoretical treatment of isolated photon cross-sections (which the collider experiments 



measure) and the related issue of the photon fragmentation function. |T^, |T^, |T8|, |T9[ This is 
relevant since the collider kinematic range is sensitive to the bremsstrahlung process, where 
a photon is radiated off an initial or final state quark. 

Analytic NLO calculations mentioned earlier [§] include the bremsstrahlung process, and 
use a leading-order (LO) fragmentation function in the region where the photon is emitted 
collinear to the quark. The implementation of the isolation cut in these calculations involves 
some approximations, and the (minor) effects of soft gluons inside the isolation cone are 
ignored. Alternatively, a Monte Carlo (MC) program for calculating NLO direct photon 
production (Baer et al. |TB|) has proved to be very useful, particularly for colliders, since 
it provides a natural way to implement isolation cuts and the soft gluon effects. The MC 
calculation yields fully-inclusive cross sections within a few % of the analytic method. For 



^The nuclear target dependence has been studied by E706 and is consistent with ° for direct photon 
productions. 
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calculations reported below, we compute the ratio of isolated-photon to total-inclusive cross 
sections using the MC method. We then use this ratio to put the collider and fixed-target 
data sets on the same footing in order to make a global comparison of data and the analytic 
calculations, and to perform new global fits for parton distributions. For the global fits, 
this procedure is also necessary since only the analytic form of the inclusive cross-section is 
efficient enough for the required repetitive calculations. 

The analytic and MC programs mentioned above both use a LO fragmentation function 
convoluted with the LO two-jet production diagrams. There is now an analytic calculation 



from Gluck et al. [19| which uses a NLO parameterization of the photon fragmentation 
function convoluted with LO two-jet production, and NLO jet production convoluted with 
the LO photon fragmentation function. In Fig. ^, we compare this calculation, as well as 
that of Baer et al. [Q, to the new CDF data. The two curves are evaluated with the same 
parton distributions (CTEQ2M) and scales (/x = pt)- The Gluck et al. calculation is 13% 
higher than Baer et al. uniformly in pf Thus the observed shape of the CDF data does not 
agree with either calculation. We note that the calculation of Ref. [|19| used GRV 0] parton 
distributions and a scale of fi = pt/2 (except for the fragmentation scale where n = (Cone- 
size) *pt = 0.7 * Pt was used). The scale change leads to a 10% normalization shift upward, 
while the different parton distributions provide a 5% shape change at photon pt = 16 GeV. 
With these parameters the total fragmentation contribution to the isolated cross section was 
estimated in Ref. [|l^ to be ~ 10%, hence any modification of this contribution is unlikely 
to explain the ~ 30% difference in shape between the CDF data and NLO QCD calculation. 
In addition, the strength of the CDF data is in their small (5%) systematic uncertainty on 
the measured shape versus pt, hence this new calculation is still insufficient to explain the 
data completely. 



4 Global Comparison of Data and NLO QCD 

The fixed-target and collider experiments described above together cover the photon xt range 
from 0.01 to 0.6. It is natural to examine all the data at once in the framework of NLO 
QCD and to utilize the combined power of this rich collection of data in a new global NLO 
QCD+parton distribution analysis. In Fig.^ we show the compilation of data sets from 
colliders and fixed target experiments compared to NLO theory (using CTEQ2M parton 
distributions and fi = Ptf^)- We see for the first time the impressive full coverage of the 
entire 0.01 < x < 0.6 range. We also see, however, that most of the data sets display a 
steeper dependence on Xt than is predicted by NLO QCD (horizontal line). We note that 
since both theory and experiment have inherent normalization uncertainties, it is difficult to 
discern an "excess" of photons at low pt from a "deficit" at high pf. However, the difference 
in slope between most data sets and the theory is quite apparent, and does not depend 
strongly on the experiment's Xt range. 

All previous global analyses used only fixed-target direct photon data in the fit. One 
obvious study is to perform a new global analysis including all the photon data in order 
to see to what extent new parton distributions, especially the gluon distribution which 
is not well constrained so far, could reduce the observed discrepancy. In this study, the 
most recent direct photon data from CDF [|15| and UA6 (both pp and pp) [ill are included 



3 



together with UA2, R806, E706, WA70 and the standard DIS and Drell-Yan data sets used 
in the CTEQ2 global analysis 0. The systematic errors for each experiment were separated 
into a normalization error (typically 10-15%) and point-to-point systematic errors that were 
added in quadrature with the statistical errors. The normalization for each experiment was 
allowed to vary from unity, constrained by the appropriate experimental errors as part of 
the minimization procedure. The theoretical scale uncertainty was taken into account by 
assigning a "theoretical error" calculated using the difference obtained with scale choices of 
and 2pt. The scale /i was then allowed to vary during the minimization process as an 
overall parameter with the associated error. 

Comparison of results from such a representative fit to the data points is shown in Fig. ^. 
The experiments are shown with their absolute normalization; in the fit they are renormalized 
downward by ~ 20%. Qualitatively, this looks similar to Fig. ^. However, the CDF data 
do come closer to the theory, while the fixed target data in particular move farther away. 
Clearly, the small statistical and systematic errors of the CDF experiment cause this data 
set to dominate the fit, hence increasing the gluon density Q) in its x-range. This 
increase is compensated by a decrease in the gluon density in the large-x region due to 
the momentum sum rule constraint .0 This demonstrates, as expected, that it is difficult to 
improve the differing slopes by new parton distributions. Furthermore, we note that it is 
also difficult to attribute these differing slopes to the photon fragmentation function since 
the ISR and fixed target data are insensitive to fragmentation. 

5 Discussion of Global Analysis Results 

What could cause an effect seen in most experiments at low pt irrespective of their x range? 
One possibility is that of "/cj smearing", since any uniform smearing on a steeply falling 
Pt distribution enhances the low pt end of the spectrum more than the high pt end. The 
question is therefore: whether the NLO QCD calculation (involving at most one non-coUinear 
initial state gluon) embodies sufficient kt broadening of the initial state partons, or whether 
additional effects due to initial state multi-gluon radiation and non-perturbative physics (i.e. 
"intrinsic kt') also need to be included? 

In order to explore this problem in a simplistic fashion, the NLO QCD prediction for the 
cross section was convoluted with gaussian functions of transverse momentum. The effect 
of such a convolution is to induce corrections to the cross section that behave as ~ cr^/pt^, 
where pt is the photon pt and a is the width of the gaussian, representing the additional 
kt in the system. The values of a needed to reproduce approximately the collider data 
sets varied from ^2-3 GeV for R806, to ^3-4 GeV for UA2 and CDF. The corrections of 
course also depend on the local slope of the NLO QCD cross sections, therefore the fixed 
target experiments are extremely sensitive to a due to their rapidly falling pt spectra. A 
width of ^1-2 GeV for E706, and slightly smaller than 1 GeV for WA70 and UA6 gave a 

^The momentum fraction carried by the gluon is constrained by the momentum sum rule to the range 
0.42-0.43 in most current global parton analyses, since the total momentum fraction carried by quarks is 
very well determined by recent precision DIS experiments. We also note that the x range of the CDF 
data (0.01-0.1) carries a large percentage of the gluon momentum, given the general shape of the gluon 
distribution. 
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reasonable overall description of all the data sets. This exercise is, of course, completely ad 
hoc: it merely serves the purpose of indicating the approximate magnitudes of the smearing 
necessary to bridge the gap between the NLO QCD calculation and experiment. The energy 
dependence of the required smearing width reveals that this possible effect is probably not 
of a simple "intrinsic" type associated with the confinement size of the hadron, but more 
likely due to multiple gluon emissions.^ 

Smearing due to multiple gluon emissions is incorporated in many QCD shower MC pro- 
grams. As another exercise, the QCD MC program PYTHIA |^ was utilized to generate the 



direct photon cross section for the CDF experiment used in our analysis. Within PYTHIA, 
there is the option to turn on/off initial state gluon radiation. We found that the ratios of 
cross sections (initial state radiation on/ initial state radiation off) show the same general 
pattern for the CDF data as seen in Figures ^ and ^. The MC program uses only LO hard 
matrix elements, and the relation of kt smearing from the multiple gluons in this calculation 
to the fixed order NLO cross-section is not clear. However, it is obvious that physical effects 
due to multiple radiation have the general characteristics of the observed cross-section. 

Theoretically, multiple gluon emission is the main physics effect treated systematically in 
the "pt resummation" formalism of vector-boson production (continuum Drell-Yan, W, and 
Z production) ||2^ . Although, in principle, pt resummation is not needed for inclusive direct 
photon production (because it is not strictly speaking a "multiple large scale" problem), it 
is intriguing to note that the resummed cross-section formula of Collins et al. [^] requires 
a "non-perturbative" broadening factor which bridges the pt physics at the confinement and 
perturbative scales, and this broadening is inherently energy dependent. It is well worth 
exploring, both theoretically and phenomenologically, any possible relations between this 
result and the question at hand. 

From an experimental point of view, one can investigate kt effects directly by examining 



the Pt imbalance of diphoton states. WA70 |23| compared the distribution of diphoton system 



Pt with that predicted by NLO QCD smeared with a gaussian kt broadening. They found 
that the typical width of 0.34 GeV intrinsic parton transverse momentum was insufficient 
to explain their data, while an additional broadening of 0.91 GeV gave a much improved 
fit. This extra broadening for diphotons is close to the amount that gave a reasonable 
overall description of the fixed target single photon measurements mentioned earlier. A 
similar situation exists for CDF diphotons ||2^. Fig. ^ shows the diphoton system pt from 



the data, NLO QCD and PYTHIA. The small statistics of the data do not provide strong 
discrimination at this time, but the important point is that the distribution from PYTHIA 
is significantly broader than NLO QCD. Once again if the NLO QCD prediction is smeared 
with a gaussian, a width of ~3 GeV is needed to reproduce PYTHIA, and this is close 
to what was needed to bring the CDF single photons in agreement with NLO QCD. The 



^In the previous analysis in ref. mentioned earlier, there was no conclusion of the need for an additional 
contribution from initial state kt- However, some elements were present in that analysis that might have 
masked this effect: (i) they emphasized the WA70 data, which are qualitatively consistent with NLO QCD 
without additional smearing; (ii) the low pt data points from the ISR (where the deviation is most apparent) 
were removed in the analysis because the scale optimization failed; and (iii) the discrepency between the 
low Pt data points of UAl & UA2 and theory was attributed to uncertainties of photon fragmentation 
contribution and isolation cuts at the time (recent progress reviewed in Sec.^has considerably reduced these 
uncertainties). 
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agreement in the magnitude of the broadening necessary for the diphotons and the single 
photons may be fortuitous, since the production processes are different. But the important 
point is that the diphotons indicate a potential problem with kt, complementing the inclusive 
photon trends. 



6 Conclusions and Future Prospects 

In conclusion we have performed a global study of direct photon experiments spanning a 
wide range of parton x values. We see a pattern of deviation between the measured slopes 
versus photon pt and the corresponding NLO QCD calculations. The deviations occur in 
several different experiments in different parton x ranges with widely varying sensitivities to 
photon fragmentation functions and theory scales. Thus it is difficult to invoke new parton 
distribution functions or a different treatment of the fragmentation process to explain all the 
data. This is verified by a new global parton distribution fit using all the direct photon data 
described in this paper. One possibility is the NLO QCD theory needs to be supplemented by 
additional kf broadening from non-perturbative physics and/or initial state gluon radiation. 
A similar phenomenon has been seen in studies of the relative pt of diphoton final states 



and this premise has been considered for some time |25]. We note that the same effect 



should also appear in low pt measurements of other processes such as dijet and heavy flavor 



production. Low pt dijet production typically has substantial systematic uncertainties |26 
but the effects of kt have been extensively studied in fixed target heavy flavor production p7[| , 
also with the conclusion for the need for additional kt broadening. On the theoretical front, 
refinement of the treatment of kt physics could take many different forms, from the analog of 
"non-perturbative functions" in Drell-Yan pt resummation formalism (cf. previous section) 
and the fc^-dependent parton distributions associated with the new "high-energy factorization 
theorem" p8|, to the interface of the parton shower calculations and NLO QCD (as performed 



for W production in p9|). 

We clearly need a better understanding of the underlying physics before low-pt photon 
measurements can be confidently relied upon as the primary source of information on gluons 
in a global QCD analysis. For the time being, the better determination of the elusive 
G{x, Q) must depend on a variety of sources of information including DIS, direct photon, 
and differential jet cross-sections. 
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Figure 1: Direct photon production data from WA70 is compared to the NLO QCD fit of 
ABFOW using optimized scales. The fractional difference is shown in order to display details 
of the comparison, (see text) 
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Figure 2: Fixed target direct photon experiments WA70, UA6, and E706 are compared to 
a NLO QCD results. Fractional differences between data points and the CTEQ2M fit are 
shown, (see text) 
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Figure 3: Direct photon data from the CDF experiment is compared to NLO QCD predic- 
tions, (see text) 
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Figure 6: Diphoton system pt distribution for PYTHIA, NLO QCD, and CDF data, (see 
text) 



14 



